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Although downbursts have been identified as the major cause of a number of aircraft takeoff and 
landing accidents, only the 1985 Dallas/Fort Worth (DFW) and the more recent (July 1994) Charlotte, 
North Carolina, landing accidents provided sufficient onboard recorded data to perform a comprehensive 
analysis of the downburst phenomenon. The first step in the present analysis was the determination of 
the downburst wind components. Once the wind components and their gradients were determined, the 
degrading effect of the wind environment on the airplane’s performance was calculated. This wind-shear- 
induced aircraft performance degradation, sometimes called the F-factor, was broken down into two 
components F, and F 2 , representing the effect of the horizontal wind gradient and the vertical wind 
velocity, respectively. In both the DFW and Charlotte cases, was found to be the dominant causal 
factor of the accident. Next, the aircraft in the two cases were mathematically modeled using the longi- 
tudinal equations of motion and the appropriate aerodynamic parameters. Based on the aircraft model 
and the determined winds, the aircraft response to the recorded pilot inputs showed good agreement with 
the onboard recordings. Finally, various landing abort strategies were studied. It was concluded that the 
most acceptable landing abort strategy from both an analytical and pilot’s standpoint was to hold constant 
nose-up pitch attitude while operating at maximum engine thrust. 


Nomenclature 


a x , a v , a z 

= body-axis accelerations 

Cl, Cp 

= lift, drag coefficients 

c La , c,„ a , c„ a 

= stability derivative 

c m 

= pitching moment coefficient 


= side-force coefficient 

D 

= drag 

F , F„ F 2 

= F-factor and its components 

h, h e 

= geometric and energy altitude 

/,, 

= moment of inertia about y axis 

4 

= pitch rate 

T 

= thrust 

U, W 

= body-axis velocity components 

V 

= true airspeed vector 

w 

= wind velocity vector 

W„ Wy, W A 

= north, east, and upward wind vector 


components 


= horizontal tailwind 

x, y , h 

= aircraft position (Earth frame) 

QLj ^vanc 

= angle of attack, vane angle 

P 

= side-slip angle 

y a 

= air mass flight-path angle 

®flap 

= flap deflection angle 
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8 h = stabilizer deflection angle 

j Jj w = wind-axis heading angle 

ip, 6, 4> ~ body-axis Euler angles 

I. Introduction 

W IND shear has been a major cause of aircraft accidents 
throughout aviation history. In the U.S., between 1964 
and 1985, wind shear has contributed to at least 26 civil ac- 
cidents involving more than 500 fatalities. 1 In the five years 
following 1986, airline pilots have reported 96 severe turbu- 
lence incidents, 2 some of them involving downbursts. 

Two common types of airplane encounters with wind shear 
are turbulence encounters at high cruising altitudes and down- 
burst penetrations during takeoff and landing. Downbursts oc- 
cur near the ground beneath storm cells and may be accom- 
panied by heavy rain (wet downbursts) or in some cases, if 
the atmosphere is very dry, may occur in the virga beneath the 
storm cell (dry downbursts). 3 Downbursts can also be classified 
as microbursts or macrobursts according to the diameter of the 
outburst winds. 3 Since the first recognized downburst accident 
in 1956, there have been at least eight other major encounters. 
A tabulation of downburst accidents/incidents is shown in Ta- 
ble 1. 

Previously, in cases of clear air turbulence incidents, data 
from the onboard Digital Flight Data Recorder (DFDR) have 
been used in conjunction with ground-based radar tracking 
(ATC radar) to extract the components of the turbulent 
winds. 4-6 Bach and Wingrove 6 provided the wind estimations 
used in the preliminary investigation of the Delta L- 1011 ac- 
cident that occurred during landing at the Dallas/Fort Worth 
airport. Fujita 7 performed DFW downburst studies. As previ- 
ously mentioned, only two of the nine downburst accidents 
listed in Table 1 involved airplanes equipped with DFDRs, i.e.. 
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Table 1 History of downburst accidents 


Location 

Date 

Airplane 

Takeoff/ 

landing 

Fatality/ 

passenger 

Kano, Nigeria 

June 24, 1956 

Boac Arogaut 

Takeoff 

32/45 

Pago Pago, Somoa 

January 30, 1974 

Pan Am 806 

Landing 

96/101 

JFK, New York 

June 24, 1975 

B-727 (Eastern 66) 

Landing 

112/124 

Denver 

August 7, 1975 

Continental 426 

Takeoff 

0/134 

Doha, Qatar 

May 14, 1976 

Royal Jordian 600 

Landing 

45/64 

Philadelphia, PA 

June 23, 1976 

Allegheny 121 

Takeoff 

0/106 

New Orleans, LA 

July 9, 1982 

B-727 (Pan Am 759) 

Takeoff 

152/161 

Dallas Fort Worth, TX 

August 2, 1985 

L-1011 

Landing 

137/163 

Charlotte, NC 

July 2, 1994 

DC-9 

Landing 

37/57 


the DFW L-1011 landing accident in 1985 and the 1994 Char- 
lotte, North Carolina DC-9 accident that occurred during an 
attempted landing abort. Thus, the lack of an adequate data- 
base restricted the present analysis to these two accidents. 

This article analyzes in detail the two previously mentioned 
encounters. Wind components in both cases were estimated. 
The degradation of the aircraft’s climb capability and ability 
to maintain airspeed, sometimes called the F-factor, was cal- 
culated and discussed. Also, flight simulations were performed 
to validate the aerodynamic model and to examine various 
landing abort strategies. 

II. Wind Estimation Methods 

Before the introduction of the flight data recorder (FDR) and 
the cockpit voice recorder (CVR), aviation accident investi- 
gators had to rely on examination of the wreckage, meteoro- 
logical data, and interviews of survivors and eyewitnesses. 
This method left a lot of questions unanswered. The early 
FDRs recorded only four channels of data: 1) indicated air- 
speed, 2) vertical acceleration, 3) pressure altitude, and 4) 
magnetic heading, along with radio transmission times. Al- 
though an improvement in the investigation of the wind-shear- 
related accidents, the information provided was so sparse that 
investigators were unable to reconstruct, with any precision, 
the wind environment surrounding the accident (e.g., Eastern 
Airlines Fl-66, John F. Kennedy Airport, New York, and Pan 
Am Fl-759, New Orleans, Louisiana). 

Digital flight data recorders were made mandatory on all 
large airliners certified after 1969. The extensive set of varia- 
bles recorded on the DFDR along with ATC radar tracking 
data now provided a sufficient data set to accurately recon- 
struct the wind environment experienced by the airplane. 4,5 

The wind vector W is defined as 


W = V, - V (1) 

where V, is the aircraft velocity in an Earth-fixed frame, and 
V is the vector of the true airspeed relative to the moving air 
mass (Fig. 1). The three components of W are given by the 
following equations: 


W x = x — V cos i p w cos y a 

(2) 

W y = y — V sin i p w cos y a 

(3) 

W h = h — V sin y a 

(4) 


where y a , the air mass flight path angle, and i p w , the wind axis 
heading angle, are calculated from the following equations 4 : 

sin y a = cos a cos (3 sin 8 — C cos 8 (5) 


tan(i/v ~ ip) = 


sin /3 cos <p — sin a cos /3 sin <p 
cos a cos (3 cos 8 + C sin 6 


( 6 ) 



Fig. 1 Velocity vector: a) plan and b) side views. 
The angle 13 is calculated using the following equation: 


(3 = (l/C Vp )(tf,C - C ya 8 r ) (8) 

where C„. is the weight coefficient, C y$ and C Vg are side-force 
derivatives, and S r is the rudder-deflection angle. 

The velocity components in an Earth-fixed coordinate sys- 
tem may be calculated, using the following equations 8 : 

x = a x cos 8 cos ip + a v (sin (p sin 8 sin ip — cos (p sin ip) 

+ a.( cos <p sin 9 cos ip + sin cp sin ip) (9) 

y = a x cos 8 sin i/> + « v (sin <p s ‘ n # sin ip + cos (p cos ip) 

+ a. (cos <p sin 9 sin ip — sin (p cos «A) (10) 


(7) 


h = a, sin 8 — (a y sin <p + a. cos <p)c os 9 — 1 


C = sin a cos [3 cos <p + sin (3 sin <p 


( 11 ) 
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These equations are integrated successively to provide the in- 
ertial velocities x , y, and h and the inertial positions x, y, and 
h. Biases and initial conditions are determined by matching the 
calculated ground positions to ATC radar tracking data and the 
calculated altitude to the DFDR recorded barometric alti- 
tude. 4- ' 6 

There are two methods of determining the angle of attack. 
In the DFW case, a direct measurement of the angle of attack 
was made with a vane. In this case, with a flap setting of 33 
deg, the true a is given by the calibration equation: 

a = 0.535a vanc + 3.72 deg (12) 

In the Charlotte case, where no direct measurement of the an- 
gle of attack was made, the angle of attack was determined 


from the lift curve (i.e., C L vs a) at the specified flap setting 
using the following equations: 


n zw = —a z cos a + a x sin a 

(13) 

Cl = n z wW/C}S = n Z wC w 

(14) 

a = (CJC La ) + a 0 

(15) 


where n zw is the wind-axis vertical load factor. In this case, 
when the pilot decided to abort the landing, the flaps of the 
DC-9 were retracted from 40 to 15 deg. In performing the 
calculation of a, since no recording of flap deflection was 
available, it was assumed that the flap retraction was started 
15 s before ground impact and that the flap deflection de- 
creased at a constant rate from 40 to 15 deg over a period of 
11s. 



Fig. 2 DFW wind (1-s interval) 





S70.0 57S.0 580.0 S8S.0 590.0 595.0 600.0 
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100 feet/sec 


Fig. 3 DFW wind (1/4-s interval). 
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100 feet/sec 

Fig. 5 Charlotte wind (1/4-s interval). 


III. Derived Winds 

Vector plots of the wind components for the DFW case are 
shown in Figs. 2 and 3 in both plan and elevation views. The 
Charlotte wind vectors are displayed in Figs. 4 and 5. As ob- 
served, during downburst encounters, the airplane, as it enters 
the downburst, first encounters strong headwinds that rapidly 
change into tailwinds. All the while, during the downburst pen- 
etration, the aircraft is subjected to downdrafts of variable in- 
tensity. 

Entrained vortices are sometimes observed in downbursts. 
They occur near the ground in the outflow along the perimeter 
of the vertical column of the downburst. The vortices, de- 
pending on their strength, can be very damaging to the air- 
craft’s performance. These entrained vortices can be seen in 
both the DFW and Charlotte cases. In the DFW case the vortex 
was sufficiently strong to cause the airplane to reach a stalling 
angle of attack as it exited the vortex. This vortex was en- 
countered about 16 s before ground impact and was possibly 


the major cause of the accident. After the vortex encounter the 
airplane never recaptured the glide slope or regained a suitable 
nose-up pitch attitude. In the Charlotte case, as shown in Figs. 
4 and 5, the vortex was much smaller and less intense and 
played a minor role in the accident. 

IV. Performance Degradation, F-Factor 

The damaging effect of a downburst is based on the fact 
that when the aircraft enters a downburst it first gains airspeed 
as a result of the increasing headwind. This may prompt the 
pilot to intuitively decrease power. A few seconds later, the 
headwind decreases as the aircraft penetrates deeper into the 
downburst. As the airplane exits the core of the downburst, it 
is being pushed downward by strong vertical winds that are 
accompanied by large increases in tailwind. This sudden loss 
of airspeed and altitude places the aircraft in a precarious po- 
sition with respect to maintaining the glide slope and a suitable 
flying speed. 




F-factor 
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To quantitatively analyze the aircraft performance capability, 
a quantity, h e , called herein the energy height and defined as 
the total energy per unit weight, was used. Thus, 

h e = h + (V 2 I2 g) (16) 

where h is a measure of the potential energy and V 2 /2g rep- 
resents the kinetic energy (V is the airspeed). The rate of the 
increase in total energy when subjected to horizontal and ver- 
tical winds, W x and W h is 9 


^-^)v 

g vj 

which, in dimensionless form becomes 

fh = (T~ D) _ _ WA 

V IV \g V) 


(T - D) 

h = i L y _ 

w 


(17) 


(18) 


where W is the airplane’s weight. 

[(V^/g) — (W h !V)] is a measure of the degradation of the 
aircraft performance capability, and is sometimes called the F- 
factor. 1 Thus, 


F = (WJV) - (WJV) (19) 

The F-factor may be broken down into components F, and F 2 , 
representing the effect of the changing headwinds and the 
downdrafts, respectively. Thus, 

F, = WJg (20) 

F, = -(WJV) (21) 

From the definition, it can be seen that the larger F becomes, 
the more restricted the airplane’s climb performance will be- 
come. When the F-factor exceeds a certain value [e.g., F > 
(T — D)/W ] over a certain finite time span, the airplane will 
be unable to maintain the energy level necessary to safely com- 
plete the landing or takeoff operation. The F-factor and its 
components for the two cases were calculated and are shown 
in Figs. 6-9. As one of the most severe downbursts ever doc- 
umented, the DFW case had an F-factor reaching a value ap- 
proaching 1.0 while traversing the core of entrained vortex. In 
the Charlotte case the downburst in terms of the F-factor is 
shown to be relatively less intense. It can be seen that, in both 





cases, Fi is significantly larger than F 2 , indicating that the hor- 
izontal wind gradient was the dominant cause of the accidents. 

V. Landing Abort Strategies 

Various landing abort strategies were examined for both the 
DFW and Charlotte accidents. A three-degree-of-freedom lon- 
gitudinal aerodynamic model was devised. The model was first 
validated by comparing the airplane’s response to the winds 
and the recorded control inputs with the DFDR recorded out- 
puts. The comparison showed good agreement. 

A. Flight Simulation 

The simulation assumed an Earth-based inertial frame and 
the aircraft as a rigid body. The linearized longitudinal equa- 
tions of the airplane are 8 


U = (1 /m)X — g sin 0 — qW 

(22) 

W- (l/m)Z + g cos 6 + qU 

(23) 

q = Mllyy 

(24) 

e = q 

(25) 
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C m = C,„ 0 + C,„ a a + C rn< {_q + q g ) + C mSf d H + C ms 8 e 

W + W x sin 0 + W h cos 0 
a tan U + W x cos 0 — W h sin 0 


(32) 

(33) 


dW„ %, 
q& " dx " V x 

q = TPV 2 

V = [(W + W x sin 0 + W h cos 0f 
+ (U + W x cos 0 - W h sin 0) 2 ] 1 ' 2 


(34) 

(35) 


(36) 


Starting with the correct initial conditions, integration of the 
previous equations forward in time using the pilot’s inputs as 
recorded on the DFDR, yields the state variables from which 
the flight path may be obtained. The simulated and recorded 
flight path and other variables such as airspeed, pitch attitude, 
and angle of attack were compared. The little discrepancy ver- 
ified the accuracy of the simulation model. 



Fig. 10 Altitude in 12-deg pitch hold landing abort (DFW). 



where U and W are ground referenced inertial velocity com- 
ponents along the body axis, and 0 and q are the pitch attitude 
and pitch rate in the Earth-fixed frame. After substituting C L , 
C D , and C,„ into the expanded form of X, Z, and M, the sim- 
ulation model for the airplane is as follows 10 : 

U = (1 /m)[(C L sin a — C D cos a)qS + X r ] — g sin 0 — qW 

(26) 

W = ( — \/m)[(C L cos a — C D sin a)qS + Z T ] + g cos 6 + qU 

(27) 

q = ( C m qSc + M r )/I„ (28) 

0 = q (29) 

where, considering the gust effect on a, q , and airspeed, V , 9 

C L = C L J_a — ar 0 ) + C L& J> Htav ^ 4- C L {q + q s ) + C Lo a (30) 

Co = C m + kC L + Co &H &H ma , (31) 



Fig. 12 Altitude in 12-deg pitch hold landing abort (Charlotte). 
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Fig. 13 IAS in 12-deg pitch hold landing abort (Charlotte). 


B. Previous Studies 

Previous studies of landing abort strategies have been made. 
Some of these studies have used potential flow to represent 
the downburst winds. 11 Such potential flow models would, for 
example, poorly represent the DFW downburst where the en- 
trained vortices, which were not included in the model, were 
a major contributor to the accident. In 1987, Bryson and Zhao 
proposed a control strategy that allowed satisfactory penetra- 
tion of a Boeing 727 through a severe downburst during take- 
off. 12 The strategy used a combination of tight climb-rate hold 
feedback and full throttle. They used with success Boeing’s 
recommendation for the B-727 of holding a constant pitch at- 
titude of 6 = 15 deg while traversing the downburst. 

C. Pitch-Hold Landing Abort 

In this article, several landing abort strategies were studied 
for the two downburst encounters. The wind inputs were those 
of the derived winds that included the entrained vortices. It 
was assumed that during the landing abort the winds were 
those affecting the accident airplane. Usually the abort altitude, 
at a given time, would be greater than that of the accident 
aircraft for which the winds were known. Consequently, since 
the horizontal wind gradients decrease with altitude, the pres- 
ent landing abort studies would, in general, be conservative 
with respect to the degrading effects of the downburst winds. 

With respect to the landing abort strategies studied, it was 
concluded that the most acceptable strategy was to apply max- 
imum thrust and hold constant pitch attitude during the abort. 
For both the DFW and Charlotte cases a pitch hold of 12 deg 
appeared to be preferable to 15 deg, the preference being based 
on the margin of airspeed above the stall. 

Figures 10 and 11 show the altitude and airspeed during a 
simulated landing abort of the L-1011. The abort was started 
23 s before impact, with the application of full engine thrust 
and a 12-deg pitch attitude hold. Figures 12 and 13 show the 
results of a landing abort of the Charlotte DC-9. The pitch 
attitude was held constant at 12 deg and throttle was set for 
maximum climb thrust. The abort was started 23 s before im- 


pact. In both abort cases satisfactory altitude and airspeed were 
maintained. 

VI. Conclusions 

The wind components along the flight path were first deter- 
mined and when plotted in vector form displayed the salient 
characteristics of a downburst. The aircraft, in both cases, upon 
exiting the downbursts, encountered vortices entrained in the 
outflow. In the DFW case the vortex was very intense and 
became a major contributor to the accident. 

To quantitatively measure the wind-induced performance 
degradation, the F-factor and its two components, F, and F 2 , 
were calculated. When plotted in component form, it was 
shown that, in both cases, Fj was much larger than F 2 . Con- 
sequently, the rapid shift from a headwind to a tailwind, as 
quantified by F,, was the major contributor to the accidents. 

Among the many simulations studied of landing aborts, it 
was concluded that the most acceptable strategy was one of 
holding constant pitch attitude while having the engines set at 
maximum climb thrust. 
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